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Abstract 



The geometric theory of Lie systems is used to establish integrability 
conditions for several systems of differential equations, in particular some 
Riccati equations and Ermakov systems. Many different integrability cri- 
teria in the literature will be analysed from this new perspective, and 
some applications in physics will be given. 

1 Introduction 

in 

Non- autonomous systems of first-order and second-order differential equations 
appear in many places in physics. For instance, Hamilton equations are systems 
of first-order differential equations, while Euler-Lagrange equations for regular 
Lagrangians are systems of second-order differential equations. A system of 
' second-order differential equations in n variables of the form x l = F % (x,x,t), 

with i = 1, . . . , n, is related with a system of first-order equations in In variables: 



x — v 

\v i = F l (x,v,t) ' i = 1 ,---n- (1) 

Therefore, it is enough to restrict ourselves to study systems of first-order dif- 
ferential equations. From the geometric viewpoint a system 

x i =X i (x,t), i=l,...,n, (2) 

is associated with the i-dependent vector field X = X z (x, t) d/dx 1 whose integral 
curves are determined by the solutions of the system. 

Unfortunately, there is no general method for solving such equations. Rele- 
vant questions about integrability are how to find a particular solution (deter- 
mined by a;(0) = xo), or a r-parameter family of solutions, or even the general 
solution (a n-paramcter family of solutions). Finally, when is it possible to find 
and how to determine a superposition rule for solutions? 
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We shall understand that to find a solution means to reduce the problem to 
carry out some quadratures. For instance, the general solution of the inhomo- 
geneous linear differential equation dx/dt — bo(t) + b\{t)x can be found with 
two quadratures and it is given by 

x{t)=ex_p^J bx(s)ds^j x ^x a + J b (t') exp f- J b x (a) dsj dt^j . 

Actually, when the systems we are dealing with are linear, there is a linear 
superposition principle allowing us to find the general solution as a linear com- 
bination of n particular solutions. For instance, for the harmonic oscillator with 
a i-dependent angular frequency u)(t): 



-uj 2 (t)x 



X = V 

v = —u) 2 (t)x 



whose solutions are the integral curves of the ^-dependent vector field X = 
vd/dx — oj 2 (t)x d/dv, if we know a particular solution, the general solution can 
be found by means of one quadrature and if we know two particular solutions, 
Xi and x 2 , the general solution is a linear combination (no quadrature is needed) 
x(i) = ki xi (t) + &2 x 2 (t) . 

There are systems whose general solution can be written as a nonlinear 
function of some particular solutions. For instance, for Riccati equation: if a 
particular solution is known, the general solution is obtained by two quadratures, 
if two particular solutions are known the problem reduces to one quadrature and, 
finally, when three particular solutions are known, x\,x 2 and X3, the general 
solution can be found from the cross-ratio relation 

X — X\ X3 — X\ , 
: = k , 

x - x 2 x 3 - x 2 

which provides us a nonlinear superposition rulepQ. 

There also exist cases in which we can superpose solutions of one system for 
finding solutions of another one. We have seen one example: the general solution 
of the inhomogeneous linear equation can be written as x(t) = x±(t) + C xo(t), 
where xo(t) is a solution of the associated homogeneous equation and Xi(t) is a 
particular solution of the inhomogeneous linear one. 

Milne-Pinney equation[2] x = —u 2 (t)x + k/x 3 is usually studied together 
with the time-dependent harmonic oscillator y + uj 2 (t)y = and the system is 
called Ermakov system. Pinney showed in a short paper [5] that the general 
solution of the first equation can be written as a nonlinear superposition of 
two solutions of the associated harmonic oscillator. All these properties can 
be better understood in the framework of Lie systems, conveniently extended 
in some cases to include systems of second-order differential equations. These 
systems have a lot of applications not only in mathematics but also in many 
different branches of classical and quantum physics. 

Let us look for systems admitting a (maybe nonlinear) superposition rule. 
The main result was given by Lie H] : 

Theorem: Given 0) a necessary and sufficient condition for the ex- 
istence of a function 4> : R Tl ( m + 1 ) — > W 1 such that the general solution is 
x = $(ai(i), . . . , x/ m y, k\, . . . , k n ), with {x^ | a — 1, ...,m} being a set of 
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particular solutions of the system and k\, . . . , fc„, are n arbitrary constants, is 
that the system can be written as 

d-r 1 

— = z\t)ti(x) + --- + z r (t)e r (x), 

where Z 1 , . . . , Z r , are r functions depending only on t and a = 1, . . . ,r, 
are functions of x = [x , . . . , x"), such that the r vector fields in 1™ given by 

n 

X a = ~y^£a( xl i ■ ■ ■ ,x n )d/dx l , a — 1, ... ,r, close on a real finite- dimensional 
Lie algebra, i.e. the X a are Li. and there are r 3 real numbers, c a f3 J , such that 

r 

[X a , Xp\ — ^^c a p 1 X 1 . The number r satisfies r < m,n. 

7=1 

The condition in the Theorem is that X{x,t) can be written as X(x,t) = 

r 

'^^Z a (t)X a (x) 1 with X a as mentioned above. 

a=l 

Non-autonomous systems corresponding to such t-depcndcnt vector fields 
will be called Lie systems. One instance of Lie system is the Riccati equation^ 

= b 2 (t) x 2 (t) + 6i(t) x{t) + b (t) , (3) 

for which m ~ 3 and the superposition principle comes from the relation 

x — X\ X3 - X\ _ ^ _ kxi(x 3 ~ x 2 ) + x 2 (xi - x 3 ) 

x - x 2 ' x 3 - x 2 k(x 3 - x 2 ) + (xi - x 3 ) 

The associated Lie algebra is generated by Xo, X\ and X 2 given by 

d d 2 d 

Xq = tt- , X\= x — , X 2 = X — — , 

ox ox ox 

which close on a s[(2,M) 3-dimensional real Lie algebra, because 

[X ,X 1 ]=X , [X ,X 2 ] = 2X 1 , [X 1 ,X 2 ]=X 2 . (4) 

The timc-dcpcndcnt harmonic oscillator is also an example of physical rele- 
vance. It is described by a Hamiltonian 

H=^^- + lm(t)^(t)x\, 
2 m(i) 2 

which gives rise to the dynamics defined by the t-dependent vector field 

vl \ 1 ° I \ 2/ \ ° 

X(x,p,t) = — — p- m(t)uj (t)x—. 

m(t) ox op 

If we consider the set of vector fields 

Y 9 Y 1 ( 9 9 \ Y 9 

ox 2 V ox op J op 
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which close on a sl(2,R) Lie algebra with the same commutation relations as 
(H]) , the corresponding i-dependent vector field X can be written as a linear com- 
bination X(-,t) — m(t)iu 2 (t) X2(-) + (l/m(t)) Xo(-), i.e. it is a linear combina- 

2 

tion with i-dependent coefhcients X(-,t) = ^^b a (t)X a (-) with bo(t) — l/m(t), 

a=0 

h(t) = and b 2 (t) = m(t)uj 2 (t). 

The prototype of Lie system is the time-dependent right-invariant vector 
fields in a Lie group G. Let {ai, . . . ,a r } denote a basis of T e G. A right- 
invariant vector field X R is one such that X R (g) = R gSfe X R (e). Define X R by 
Xaid) — Rg*e&a- The independent right-invariant vector field 

r 

X(g,t) = -Y,b a (t)X*(g) , 

defines a Lie system in G whose integral curves are solutions of the system 
g = —J2a=iba(t)X R (g), and when applying R g -i to both sides we see that 
g{t) satisfies 

r 

Rg-l(t)*g(t)9(t) = 6 "W a « € T eG . (5) 

a=l 

Let H be a closed subgroup of G and consider the homogeneous space M — 
G/H. Then, G can be seen as a principal bundle over G/H: (G,t,G/H). 
The X R are r-projectable on the corresponding fundamental vector fields of 
the left-action A : (g,g'H) eGxM-t (gg'H) G M given by -X a = -X Ua 
with T*gX R (g) = —X a (gH), the projected vector field in M will be X(x,t) = 
J2a=i^a(t)X a (x), and its integral curves are the solutions of the system of 
differential equations: x — ~Y^ r a= ib a {t)X a {x). The solution of this last system 
starting from x$ is x{t) = X(g(t), xq), with g(t) being the solution of ((SJ) such that 
g(0) = e. This means that solving such Lie system in G we are simultaneously 
solving the corresponding problems in all its homogeneous spaces. 



2 SODE Lie systems 

A system of second order differential equations can be studied through the corre- 
sponding system of first-order differential equations as indicated in (fT]) . We call 
SODE Lie systems those for which the associated first-order one is a Lie system, 
i.e. it can be written as a linear combination with t-dependent coefficients of 
vector fields closing on a finite-dimensional real Lie algebra. An example is the 
1-dimensional harmonic oscillator with time-dependent frequency, but the same 
is true for the 2-dimensional isotropic harmonic oscillator with time-dependent 
frequency, with an associated vector field 

d 2 d d 2 1 \ ^ 

A = «i- u [t)xi ^ h w (t)x 2 , 

ax\ av\ ox 2 0V2 

which is a linear combination X = X2 — oj 2 {t)X\ with 

d d v d d 

Xl=X X - hl27T . X 2 =Vl- hl>27 , 

ovi ov 2 ax\ 0x2 
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and then they close once again on a Lie algebra isomorphic to s[(2,R): 



[X 1 ,X 2 ]=2X 3 , [X 1 ,X 3 ]=-Xi, [X 2 ,X z ]=X 2l (6) 

with 

X --( — - — — - — 

3 2 \ 1 dx\ 1 dv\ 2 dx2 2 dv2 

The search for the superposition rule for a Lie system consists on looking for 
enough number of first integrals, independent of the time-dependent coefficients, 
in an extended space in which we consider several replicas of the given vector 
field. 

The 2-dimensional case admits an invariant F given by the first integral 
F(xi, Xs, Vi, v 2 ) — x\V2 — X2V1, which can be seen as a partial superposition 
rule[5]- Actually, if x\(t) is a solution of the first equation, then we obtain 
for each real number k the first-order differential equation for the variable x 2 , 
xi(t) d,X2/dt = k + Xi(t)x2, from where X2 can be found to be given by x 2 (t) = 
k'xi(t) + kxi(t) f x^ 2 (£) d£. With three copies of the same harmonic oscillator, 
i.e. Xi and X 2 given by 

d d d d d d 

Xi=vi- \-v 2 - \-v—, X 2 = xi- hx 2 - ha;—, 

oxi 0x2 ox av\ 0V2 ov 

there exist two independent first integrals F\(x\,X2,x,V\,V2,v) = xv\ — x\v 
and F2(xi,X2,x,vi,v 2 ,v) — XV2 — X2V, from where we obtain the expected 
superposition rule: 

x = ki x\ + k 2 x 2 , v — fci Vi + k 2 v 2 ■ 

Another interesting non-linear example is the Pinney equation, the second 
order non-linear differential equation: 

2 / \ ^ 

X = -Id (t)x H — T , 
X 6 

where A; is a constant, with associated i-dependent vector field 

d ( 2 / \ k \ d 
X = V dx + {-" {t)X+ xt)d- V > 

which is a Lie system because it can be written as X = L2 — uo 2 {t)Li : where 
L 1 := xd/dv and L2 = (k / x 3 ) d / dv + v d / dx generate a three-dimensional real 
Lie algebra isomorphic to sl(2,R) with nonzero defining relations similar to © 
with L 3 = (1/2) (xd/dx-vd/dv). 

Note that this isotonic oscillator shares with the harmonic one the property 
of having a period independent of the energy, i.e. they are isochronous, and in 
the quantum case they have a equispaced spectrum. The fact that they have 
the same associated Lie algebra means that they can be solved simultaneously 
in the group SX(2,R) by the same equation 

Rg-^*g9 = u 2 {t) a% - a 2 , g(0)=e. 



5 



3 Ermakov systems 



We can consider the generalised Ermakov system given by: 

1 



1 



■f(y/x)-u>*(t)x 
g(y/x)-ui 2 (t)y 



y = ~3 



which when f(u) — k and g(u) = reduces to the Ermakov system. 
This system is described by the ^-dependent vector field 



d 



d 



dx dv v 



1 



d 



dv x 



1 



V 



d 



dv v 



which can be written as a linear combination X = N 2 — u 2 (t) Ni , where Ni and 
N 2 are the vector fields 

that generate a three-dimensional real Lie algebra isomorphic to sl(2,R) with 
a third generator 



d 



d 



d 



d 



N z = 77 x 7> v ^Ti 1" yj, v v^~ 

ox dv x dy dv y 



There exists a first integral for the motion, F : M 4 —> K, for any o-i 2 (t), which 
satisfies NiF = for i = 1, . . . , 3, but as [Ni, N 2 ] — it is enough to impose 
NiF = N 2 F = 0. The condition NiF = 0, implies that there exists a function 
F : M. 3 — > K such that F(x, y, v x ,v y ) = F(x, y, £ = xv y — yv x ). Then using the 
method of the the characteristics in condition N 2 F = 0, we can obtain the first 
integral: 



1 [ x l y \ 1 / 1 \ / 

F(x,y,v x ,v y ) = — (xVy - yv x f + / -— / ( - J +ug ( 



du . 



For the Ermakov system with f(l/u) = k and g(l/u) = we obtain the known 
Ermakov invariant 

F(x,y,v x ,v y ) = - {-) +-(xv y -yv x ) 2 

We can now consider a system made up by a Pinney equation with two associ- 
ated harmonic oscillator equations, with associated ^-dependent vector field 



d 



d 



d k d 



X = v x — + v y — +v z — + -r- lu (t) x- h y-^ h z— 



1 dy z dz y 3 dv. 



d 



d 



d 



dv x dv v 



dv. 



which can be expressed as X = N 2 — oj 2 (t)Ni where Ni and N 2 are: 



Ni=y 



_d_ 

dv v 



d d did d 

x- hz — , N 2 =v y — + 

dv z 



dv x 



dy y 3 dv y x dx 



+ v z 



dz' 
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These vector fields generate a three-dimensional real Lie algebra isomorphic to 
s[(2,R) with the vector field iV 3 given by 



, ■ d d d d 

^3=o x a v x7> >" Vji v vTi !" z 7> v = a 

ox Ov x oy ov y oz ov z 



In this case there exist three first integrals for the distribution generated 
by these fundamental vector fields: The Ermakov invariant I\ of the subsystem 
involving variables x and y, the Ermakov invariant I2 of the subsystem involving 
variables y and z, and finally, the Wronskian W of the subsystem involving 
variables x and z. They are given by W = xv z — zv x , 

h = 7; ( {V v x -xv v ) 2 +c ( - ) ], I 2 = \{^V v z ~ zvyf + c (J 

In terms of these three integrals we can obtain an explicit expression of y in 
terms of x, z and the integrals I\, 7a, W: 




V= W i 12 ^ +hz2± V4/i/2 - cW 2 xz 



1/2 



This can be interpreted as saying that there is a superposition rule allowing 
us to express the general solution of the Pinncy equation in terms of two inde- 
pendent solutions of the corresponding harmonic oscillator with time-dependent 
frequency. 



4 The reduction method and integrability crite- 
ria 

Given an equation (0) on a Lie group, it may happen that the only non- vanishing 
coefficients are those corresponding to a subalgebra f) of Q and then the equa- 
tion reduces to a simpler equation on a subgroup, involving less coordinates. 
An important result is that if we know a particular solution of the problem 
associated in a homogeneous space, the original solution reduces to one on the 
isotopy subgroup. 

One can show that there is an action of the group Q of curves in G on the 
set of right-invariant Lie systems in G (see e.g. [5] for a geometric justification), 
and we can take advantage of such an action for transforming a given Lie system 
into another simpler one. 

So, if g(t) is a solution of the given Lie system and we choose a curve g'(t) 
in the group G, and define a curve g(t) by ~g(t) — g'{t)g(t), then the new curve 
in G, g(t), determines a new Lie system. Indeed, 

r 

Rg(t)-i*g(t)(g(t)) = Rg'-i( t )*g<(t)(g r {t)) - ^2 b a (t)Ad (g'{t))& a , 

a=l 

which is similar to the original one, with a different right-hand side. Therefore, 
the aim is to choose the curve g'(t) in such a way that the new equation be 
simpler. For instance, we can choose a subgroup H and look for a choice of g'(t) 
such that the right hand side lies in T e H, and hence ~g(t) e H for all t. This can 
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be done when we know a solution of the associated Lie system in G/H allows 
us to reduce the problem to one in the subgroup i?[7]- 

Theorem: Each solution of 0) on the group G can be written in the form 
g(t) = gi(t)h(t), where gi(t) is a curve on G projecting onto a solution gi(t) 
for the left action A of G on the homogeneous space G/H and h(t) is a solution 
of an equation but for the subgroup H , given explicitly by 

(R h - Uh h)(t) = -Ad^r'W) + (R gT ^ gi gi)(t)J e T e H . 

This fact is very important because one can show that Lie systems associated 
with solvable Lie algebras are solvable by quadratures and therefore, given a Lie 
system with an arbitrary G having a solvable subgroup, we should look for a 
possible transformation from the original system to one which reduces to the 
subalgebra and therefore integrable by quadratures. 

By the last Theorem there always exists a curve in G that transforms the 
initial Lie system into a new one related with solvable a Lie subgroup of G. 
Nevertheless, it can be difficult to find out a solution of the equation in M that 
determines this transformation. Then, to be able to obtain one is more inter- 
esting to suppose also that this transformation is a curve in a certain subset of 
G, i.e. a one-dimensional subgrop. It would be easier to obtain a transforma- 
tion but it may be that such a transfrormation does not exist. In summary: 
The conditions for the existence of such a transformation of a certain form are 
integrability conditions for the system. 

We could choose for showing this assertion a particular example: Riccati 
equation. One can find in the literature a lot of integrability criteria for Ric- 
cati equation [H EH HUJ , all of them particular examples of the above method 
[TT) . We can also consider other equivalent examples as the Pinney equation, 
the (generalized) Ermakov system or more relevant examples in Physics, for 
instance, time-dependent harmonic oscillators. The results obtained for one 
system are valid for the other; they are essentially conditions for the equation 
in the group, and all are examples of Lie systems associated with the same Lie 
group: SL(2, R). Consider, for instance, the Riccati equation Q. The group 
SX(2,R) contains the affine group (either the one generated by Xo and X\ or 
the one generated by X\ and -^2), which is SOLVABLE. Therefore, a transfor- 
mation from the given equation to one of this subgroup allows us to express 
the general solution in terms of quadratures. This happens when we know a 
particular solution x\ of the given equation: x = x\ + z, what corresponds to 
choose 

««)-(! T ) 

reduces the equation to dz/dt = (2 62 x\ + b\)z + b-2 z 2 . The reduction by the 
knowledge of two or three quadratures has also been studied from this perspec- 
tive and similarly for the Strelchenya criterion [TTJj . 

Each Riccati equation can be considered as a curve in R 3 and we can trans- 
form every function in R, x(t), under an element of the group Q of smooth 



8 



SX (2, Revalued curves Map(R, SX(2,R)), as follows: 

8(A,oo) = a{t)h(t), @{A,-5{t)h{t)) = oo , 

— it *SH- 

The image x'(t) — <d(A(t), x(t)) of a curve x(t) solution of the given Riccati 
equation satisfies a new Riccati equation with the coefficients b' 2 , fo' l5 b' : 

b' 2 = 8 2 b 2 - 87 bi + 7 2 b Q + j8 - £7 , 

b[ = -2/35 6 2 + (aS + /?7)&i - 2a-fb + 8a - aS + /3j - j[3 , 
K = P 2 b 2 - a/3 bi + a 2 b + a/3 - /3a . 

This expression defines an affine action of the group Q on the set of Riccati 
equations or analogous Lie systems. 

Lie systems in SX(2,R) defined by a constant curve, a(i) = J2 a =o c a & a, are 
integrable and the same happens for curves of the form a(i) = D(t) (j2^ =0 c Q a Q ^ , 
where D is an arbitrary function, because a time reparametrisation reduces the 
problem to the previous one, i.e. the system is essentially a Lie system on a 
one-dimensional Lie group. 

We can prove the following theorem which is valid for both Riccati equation 
and any other Lie system with Lie group 5L(2,R): 

Theorem: The necessary and sufficient conditions for the existence of a 

transformation: y' = G(t)y, i.e. A(t) = ( a ^ -Iaa )> re ^ a ti n 9 the Riccati 



or^{t) 

equation f3J) with (for 6062 ^ 0, with an integrable one given by 

^=D(t)(co+c iy ' + c 2 y' 2 ), (7) 
where Ci are real numbers, Ci S R, is that c^c-i ^ 0, and: 

The unique transformation is then y' = (b2(t)co) 1 ^ 2 (bo(t)c2) ~ 1 / 2 y. 
As a consequence, given ([7]) if there are constants K, L such that 




b (t) 
b (t) 



K 



then there exists a time-dependent linear change of variables transforming the 
given equation into the solvable Riccati equation ([7]) with c\ = K,cqc 2 = L and 
D(t) is given as above. 

The existence of such constant K can be considered a sufficient condition for 
integrability of the given Riccati equation or the corresponding Milnc-Pinncy 
equation. 
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